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I 



INTeODUCTION 



Simulaticn to provide answers to difficult prcbleas is 
not new. An attempt to simulate the distribution of 
navigational errors wnen using an existing cr proposed map 
product is new. A field experiment conducted by the U. S. 
Army Combat Eevelopiaents Experimentation Command (CDEC) , 
MAFEEO III, during the period April-May 1977, was the 
motivation for a simulation approach to the distribution of 
navigational errors achieved during this experiment. This 
thesis is organized as follows to explain the simulation 
approach taken and the results of that simulation. 

A. SECTION II - 3ACKGSCUND . ' ’ 

This section is important for the reader who is 
unfamiliar with the rJAPPEO III experiment. It provides an 
overview cf the experiment tc include a discription of tne 
experiments' location, players used, the map products used, 
design and cbjectives. The results of P.APPRO III and 
analysis tectnigue CDEC utilized was provided as a departure 
point for a continuing analysis. 

E. SECTION III - DESCBIPTION OF DATA 



This section describes the data used fcr analysis. A 
description cf the raw data, generation of the data points, 
transformations used, accuracy of the data, fora fcr 
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analysis and quantity cf the data is picvided. A sunaiaiy of 
the radial navigation errors is provided as a link to the 
goal of sinulating these errors with an appropriate 
distribution . 



C. SECIICN IV - NATURE CF THE EEOBLEH 

The nature cf the protleai is a simple statenent cf the 
gcals of the thesis in relation to what was a-ccomplis bed by 
CCEC and what could be achieved by ap f roximatinc the 
distribution of the radial navigation errors. 

C. SECIICN V - TESTS CN THE CATA 

A description cf the tests performed is given tc 
determine pertinent facts about the' data. "Sirce a 
sioiulaticn of the radial navigation errors could be 
accomplished by simulation cf the components of the radial 
error, tests were performed cn the components to determina 
if these ccmponents were correlated. Since all cf the 
components could not be definitely defined as uncorrelated, 
a component approach to simulation was eliminated from 
consideration and the rerainder of the tests vere performed 
to test the homogeneity of the data points of tht radial 
errors. 



E. SECTICN VI - SIMULATED DISTRIBUTIONS 



Two distributions vere simulated tc provide a close 
approximaticn to the distribution of the radial errors. The 
Exponential simulation was performed as a result of CDEC*3 
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I 




II 



EACKGfiOUNC 



A. lEST CF KEW ANE IMfROVED HAPS AND MAP PfiCCUCTS 



1 . General 

Curing the period April-May 1977, Ibe United States 
Army Combat Developments Experimentation Command (CEEC) , 
Fort Ore, California conducted a Test of New and Imfrcved 
Maps and Map Products (Short title MAPPEO III) in order to 
evaluate the four map products used by armor personnel in 

(i) describing specific Military Graphic Information (KGI) , 

(ii) navigation of specified day and night routes to 
preselected checkpoints and 

(iii) locating specified targets during da j light and nignt 
time envirenments. 

2 • IJS^r imen t D escripti cn 

MAPPEO III was a three-part experiment. Part 1 
evaluated the ability of selected armcr personnel to 
identify MGI. Part 2 evaluated the effects on route 
planning, navigational accuracy and target location during 
daylight hours. Part 3 was a duplication of Part 2 except 
that it was conducted during periods of light levels less 
than O.OCC25 foot candles. 
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3 . 



All trials and tests were conducted at 5ort Kunter 
Liggett, California in the upper Gatilan Valle;y in 

f lat-to-rclling and rolling-to-rough terrain. 

Partic ip ation 

Thirty-six players participated in Part 1. Twelve 
Armor officers (2nd Lt) and thirteen armor mid-grade 
noncommissioned officers were available for participation in 
parts 2 and 3. Hoiever, only eight officers and eight 
noncommissioned officers actually participated in the trials 
of parts 2 and 3. A noncommissioned officer performed a 
dual role as vehicle driver and ccntr'cller during all 
trials. 

5 • Pr oduc ts 

four 1:50,0C0 nap products were involved in the 

tests : 



(i) Standard Line (Map A) . Army standard issue rap. 



(ii) Enriched Line (Map 1 ) . 
standard line map as a base, 
geographic information is added and 
shaded relief. 



A map variety using the 
Militarily significant 
utility enhanced through 
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(iii) Orthophotoiiap (Map C) . A fhoto-basec map 
proiiuct produced froE a color losaic of aerial photcciapns 
and printed using a three-color process. 

(iv) Black and White Photo (Map D) . A map prccuced 
froit a single photograph of a mosaic of photographs and 
containing a military grid system and marginal information. 



6 • Ixesrimental D esi gn 



personae 
products 
standard 
Labcrato 
to evalu 
location 
lanes in 
Along e 
end of e 
used to 
each had 
had thr 
and targ 



Part 1 of the experiment was designed to test player 
1 in thexr ability to identify MGI on the map 
. This was accomplished through the use of several 
examinations prepared by the Engineer Topographic 
ry. Parts 2 and 3 of the experiment were designed 
ate route planning, navigational accuracy and target 
s, accomplished by estanlishing fcur navigational 
the Gabilan Valley for both day and night trials, 
ach lane were six preselected checkpoints and at the 
ach lane, individual observation points (OP's) were 
locate targets. OP's at the end of lanes 1, 2 and 4 
two targets to locate. The OP' at the end of lane 3 
ee targets to locate. The lanes, checkpoints, CP's 
ets were different for day and night trials. 



7 • Obj ecti ves 

The four principle objectives of the experiment were 
to obtain data on 

(i) Kcw clearly each map product depicts types of 
MGI during daylight conditiors; 
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(ii) The usefulness of each map product in 
visualizing various terrain features and in ansi*ering 
relief-related guestioniS; 

(iii) of four different map products on route 
planning, navigational accuracy and target location curing 
daylight hours; 



(iv) The effects cf four different map products on 
navigational accuracy and target location during periods of 
low light visability (0.00025 foot candles). 



E. HESUITS Cr MAPtfC III 



1 Graph ic Information ■ 

The players answered written guesticns about MGI 
features most accurately with Wap A and least accurately 
with Map D. 



2* Cay navigation 

The Analysis of Variance (ANOVA) cross-products 
considering only the map/roura combinati-ons were 
statistically significant. A Chi-Sguare goodness of fit 
test of the sample distributions of day navigational radial 
error indicated that the closest fit was obtained by an 
exponential distribution. 
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ir/i 



3 • ta rget L ocati on 

The largest radial error coaponerts of ranee and 
deflection occurred with dap D. There was little difference 
between these error coipcnents with the other maps. 



2§tget Location 



There was no statistically significant difference 
among map prcducts. The dispersion of radial error was 
slightly larger with dap D. 

5- Niaht Ijayigaticn 

Cf the ANOVA cross-products, the trap by route and 
set of players by route were statistically significant. 
Navigation was least accurate wdth Map D and Map C cn the 
most difficult route (lane 4) . 

C. ANALYSIS TECHNIQUE 

A three-way ANCVA was performed by CDEC or the 
navigation radial error. The ANOVA considered as possible 
sources of variation the map product, the routes ard the 
sets of players. A set of players consisted cf two officers 
and two noncommissioned officers who used ail four map 
products on all four lanes. The players were assignea sets 
prior to record trials in a manner designed to minimize the 
dixfarences between sets. This assignment was accccplished 
by reviewing players* scores on the written tests and 
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evaluating their field performance during exploratory 
trials. The best players were assigned to work with the 
worst players, and the average players were grouped 
together. Eationale for establishing sets was that this 
permitted normal ANOVA techniques which requite entries in 
all data cells, the raw data for the ANOVA's were converted 
by a leg 10 transformation, which tends tc normalize and 
stabilize the variance. One of the assumptions for ths type 
of analysis is ncrmality of data. 
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III. DISCRETION Of DATA 



A. SOUfiCf 

The raw data were furnished by the United States Army 
Combat Developments Experimentation Command, Fort Crd, 
California. The data points were generated during the 
conduct of the Test cf New and Improved Maps and Map 
Products (Short title MAPPRO) Phase III, during the period 3 
through 23 May 1977 at Fort Hunter Liggett, California. 



E, NATURE Of THE DATA 



1 . 

The raw data were in the form of ten-digit 
coordinates giving iccation to the nearest meter, cedes 
identifying the lane, day/night trial, player, checkpeinr 
and map prodret used. 



■2 • Genera ti on of Data Feints 



The 

procedures, 
to each cf 
player with 



data points were generated by daylight and night 
The daylight procedures were tc assign a player 
the four navigation lanes available, present the 
a map product with six designated checkpoints 
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for a flayer to .plan his own navigation route to the 
identified checkpoints. For safety reasons during night 
trials, the players were provided the checkpoints and a 
designated route on the map product. As a player navigated 
his planned or designated route, his position was constantly 
monitored by electronic location designation devices. Cnee a 
player decided that he was at the designated checkpoint, he 
informed an his controller who activated a switch oh his 
transmitting unit which accompanied the player and 
controller in their vehicle. The act of turning on the 
transmitter designated the location of the vehicle. 

^ • Transformation of the Da t§ 

The raw data points in the form of coordinates were 
then compared to the coordinates of the designated location. 
The actual location was thea subtracted from the player* s 
location which resulted in positive and negative values on 
an x,y axis with the designated location being (0,0) on the 
axis. Positive values are errors in the north and east 
directions. Negative values are errors in the south and 
west directions. 



4 . Accuracy 

The accuracy of the radial errors calculated was determined 
by the measuring devices used. Because these calculations 
are to rhe nearest meter, the errors are probably continuous 
and identical values are due to rounding, even though the 
underlying distribution is centinuous. 

for player safety, the actual location of the player 
was continuously monitored. When a player deviated more 
than 500 meters during daylight trials, and more than 50 
meters during night trials from his selected or prescribed 
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xcute, the player was guasticned to deteraiixe if he knew his 
location. If he did net know his location, he was declared 
lest and returned to his previous checkpoint. This 
procedure, although assuring safety on the courses and being 
an cperaticnal necessity, limited the radial errors which 
would otherwise have been observed. The nuirter of everts in 
which a player was lest is an important parameter with an 
impact upon the performance of the four map products. The 
number of lost player events, however, was small in 
comparison tc the number of navigarien points used ir the 
sample. Table I gives the number of lost player events for 
each map and the total lost player events. There were 96 
checkpoints for each of the map products and 384 checkpoints 
total for each of the daylight and night trials. No record 
was made if a player recorded multiple lost player evenrs 
when trying to reach a single checkpoint. Table I lists the 
number of lost player events recorded for each map during 
the .two trials. 

TABLE I - LOST FLAYER EVENTS 



DAYLIGHT TRIALS 



MA? A 
5 



MAP B 
4 



MAP C MAP B 

5 4 

NIGHT TRIALS 



TOTAI 

18 



MAP A MA? B MAP C 



MAP C TOTAI 



1 



13 



The sample sizes for lest players 
statistical inference in the data. 
CEDC that there is no difference in 
during Eaylight Trials and that 



is toe small tc make 
The conclusions drawn by 
the effect of the maps 
maps A and B performed 
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t€tt€X aith fewer lost player events during right trials was 
subjective and probably accurate. Without the benefit of 
knowing which radial errors were achieved after a lost 
player event and the possiblility of purging these data 
points, nothing more can be said about these events and 
their effect was not considered in the remainder cf the 
analysis . 



iSI A D^ly sis 

The ten-digit coordinate of a player's location was 
split into two parts: the first five figures being 
East/West and the second five figures being North/South. 
These figures translate naturally into the x and y axis. 
Absolute miss distances from the actual checkpoints was then 
calculated by subtracting the Kcrth/i>outh. and East/Hest true 
checkpoint figures from the player's Nerth/South and 
Easr/West figures. This provided X and X readings which 
were individually sguared and then the X-sguare was summed 
with the associated Y-sguare. A sguare-rcct was then taken 
of this sum and the result was the absolute miss distance 
from the actual checkpoint, subsequently referenced as the 
Eadial Error and the data form used for analysis. 

6 . fiiia nti_^y o t Data 

There were four map products tc evaluate, four 
navigational lanes and six checkpoints on each lane. Each 
map product was used on each navigation lane resulting in 96 
readings for each map during daylight trials and 96 readings 
for each map during night trials. The total number of 
readings was 384 for daylight trials and 384 for night 
trials. 
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C. SUWaABY CF THE SADIAL EEFCfiS 
Sutimaries of the calcirlated caraineters of the Radial Errors 



are given in Table 


£ II and 


III. All 


values. 


except for 


Skewness and Kurtosis, are in esters. The 


Skewness is a 


ireasure of the 


iiv 

sy cetry 


about the 


fflea n . 


The Kurtesis 


indicates the shape of 


the density 


f* function with large 


values indicating 


sharp peak and snail 


value £ 


indicating a 


flat peak. 

TABLE il - SUMMARY CF RADIAL ERRORS 


- EAYIIGHT TRIALS 




MAP A 


MAP B 


MAP C 


MAP D 


Mean 


42.14 


44.24 


47.23 


47. 17 


Standard Deviation 


63.60 


65.18 


61 .C6 


60.24 


Skewness 


2.516 


2.687 


2.. 062 


1.831 


k urtosis 


5 .90 0 


7.802 


3.654 


2.963 


10th Quantile 


3.606 


3.606 


5.365 


3.162 


25th Quantile 


8.395 


7.176 


9.641 


7.839 


50th Quantile 


19.35 


17.73 


22.34 


17.46 


75th Quantile 


38,97 


54. 15 


51.92 


72- Cl 


SOth Quantile 


115. 8 


107.8 


138.7 


126.6 


Kaximuoi 


3 19.7 


357.8 


274.6 


282 .0 



21 



1ABL£ III 



SUMMARY CF RACIAL ERRORS 



RIGHT TRIALS 





MAP A 


MAP B 


MAP C 


MAP D 


M€aD 


46.32 


63.06 


54.60 


63.00 


Standard Deviation 


60.20 


69.76 


74.75 


60. 19 


Skewness 


2.26 8 


2.773 


3.816 


1.490 


Kurtosis 


5.28 3 


1 1.33 


20.63 


2. 253 


10th Quantile 


3. 162 


5 .385 


5.365 


4.243 


25th Quantile 


7.567 


12.67 


10.75 


14.49 


5Cth Quantile 


19.26 


45.79 


33.56 


50. 26 


75th Quantile 


61.48 


87.23 


6 2.60 


86.79 


90th Quantile 


1C9. 0 


137.4 


135.1 


161 .0 


Maximum 


299. 0 


464.1 


558.4 


505 . 1 
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IV. 



BAT UBE Cr THE BROELEjl 



CDEC utilized a transformation of the racial errors in 
order to perforai ANOVA on the data to detercine differences 
in xhe map products. The only attempt to identify the 
distribution of radial errors was a Chi-Sguare Goodness of 
Fit Test which accepted the hypothesis that the errors were 
Exponentially distributed. Since CCEC could not rack the 
effectiveness of the map products using AhCVA techniques, 
knowledge of the distribution of radial errors could provide 
fetter methods to arrive at a decision as to which map 
product performed best. The problem, therefore, is tc 
detiirmine the actual distribution of the radial errors. 
This Thesis is . a simulation approach to finding a 
distribution which best fits the actual distribution c’f tne 
navigation radial errors. Two distributions. Exponential 
and Gamma, were simulated and compared to the actual data, 
cnce a distribution for the radial errors has been 
approximated, decisions about the perfcrmance of each map 
product can te made. 
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V. 



TESTS ON THE DATA 



A. PARAMETBIC TESTS 

One afprcach tc sinulating the radial errors would be to 
simulate the X and Y errors and then convert these valtes tc 
radial errors. In crier to simulate the 1 and Y values, 
knowledge of the correlation between the X and Y values was 
required. If little of no correlation existed between these 
values and the distriouticn of these values could be 
approximated by the Normal distribution, the problem of 
simulation would be -reduced. An analysis of the X and Y 
values using a histogram was not. encouraging . in that the 
values did not appear Normally distributed. The following 
test assumes that the X and Y values are Ncmally 
distributed. The HIKED Statistical Package £Saf. 6] was 
used in the following test. 

^ Plotting 

The bivariate plotting of the X, Y data was 
performed tc observe the relative locations of the (X,Y) 
errors ir, relation to a common checkpoint and to determine 
if trends could be observed. These plots were prepared for 
each map for daylight and night trials. Mith each plot, the 
EIKED Statistical Package provided the correlation, the 
means and standard deviations of X and Y, a linear 
regression line and the residual mean squares. The results 
of this analysis are described in the fcllcwing sections. 
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2 . Correlation Bet ween 2^ and Y 

Since the ElilEE bivariate plotting cf the X anc Y is 
tased cn the assumption of Normality, tests or the 
Coefficient cf Correlation were based on the hypothesis that 
the X and Y values were independent if the absolute value of 
the Coefficient of Correlation was not too large. A 
Student's "t" test was used for the hypothesis that the 
Coefficient cf Correlation equaled zero. Ihe "t" statistic 
was calculated by: 

Let: r = Sample Coefficient cf Correlation 

N = 96 

Then: t = r (SQ3T (N - 2))/SQKT(1 - r*r) 

Ihe significance of each Coefficient cf Correlation 

determines whether in fact the X and Y .values afe 
independent (for Normally distr-ibuted X and Y values) . The 
results fcr daylight navigation trials indicated a negative 
coefficient cf correlation cf 0.282 (Significant at Alpha = 
0.003) for Hap A, a positive coefficient cf correlation of 
0.154 (Significant at Alpha = 0.07) for Map B, a negative 
coefficient cf correlation of 0.223 (Significant at Alpha = 
0.008) fcr Hap C and a positive coefficient of correlaticn 
cf 0.085 (Significant at Alpha = 0.20) for Map D. Map D 
accepts the hypothesis for any Alpha less than 0.20. Map B 
accepts the hypothesis for any Alpha less than 0.07. Maps A 
and C reject the hypothesis for reasonable Alpha values. 
This indicates that tne values fcr at least two of the maps 
are correlated under the assumption of Normality. Ihe 
results for night navigation trials all reject the 
hypothesis. 
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Otserved Trends in the Plots 

For the daylight trials, most of the (X,Y) errors 
were to the Hast and South of the actual checkpoint. Ihis 
result is net surprising, given the human tendency to drift 
to the right when attempting to go in a straight line, and 
most of the routes to the checkpoints followed a general 
asmuthe of East to Southeast. During the niebt trials, the 
same trend in West and South errors was ctserved. Eaps A 
and C shew trends in that the radial errors are in a 250 
meter card frem Northwest to Southeast. Maf E and D do not 
exhibit this trend. 

E. NONPAEAMEISIC COBREIATIOS TESTS ON X AND Y DATA 

To further investigate , the poss-ibilitj of correlation 
between the X and Y values, several Nonpar ametric Tests were 
performed. The bivariate Plotting and resultant correlation 
provided previously, assume that the X and Y values are 
Normally distributed. Since this fact has not been shewn, 
Ncnparametric tests for correlation were performed. Ail 
tests shewn below were performed by the EIMED Statistical 
Package. 
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Spgaraan Rank Co rre la ti cn ^€fji£is£t iSgearaar Fho). 





The 


results 


of 


the 


test 


were 


as follows 


1 








MAP A 


MAP B 


MAP C 


MAE C 


Daylight Tri 


als 


-0. 


275 


-0 


.015 


-0.C6S 


0.0C7 


Night 


Trials 




-0. 


209 


-C 


.324 


-0. 410 


-0.272 


Since 


the 


number 


of 


samples 


is 


greater than 30, the 


distribution 


of the 


Spe 


aroan Rho 


is 


adequately 


approximated. 


under 


the hypothesi 


3 cf 


no 


trend 


/ by 


a normal 


distribution 



with mean 2 ero and variance (V) equal tc the recipical c£ 
the sample size minus cne (V = 1/(96 - 1)). The value 
C. 01053 was then used as the variance, so the Standard 
Ceviation was 0.1026. The hypothesis cf nc trend could then 
te rejected with an Alpha = 0.10, if the absolute value of 
Spearman Bho was greater than. 0.13150. By this criteria Map 
A, during daylight trials, and all maps during the night 
trials, show correlaticn between the X and Y values. 



2* Cth^ Correlation T^sts 

The Kendall Rank Correlation Coefficient Test was 
alsc perforned and supported the findings of the Spearman 
Rho and is net included. Since a correlaticn is indicated 
between the Y and Y values, a sinpie simulaticn of the X and 
Y values was not possible and the approach was eliminated 
from consideration. 
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C. NONP AfiAtlEISIC TESTS POP. HOMCGEWEITY OF BACI4L EREOfS 

Sinc€ th€ X and Y values were not uncorrelated for all 
maps, the simulation of the radial errors ty simulating X 
and Y was discarded and emphasis was placed cn the radial 
errors. Tests were performed cn the radial error values ro 
determine if the underlying distributions fcr the radial 
errors from all map products were the same. Two 
Ncnparametric tests were used; The Kolmogcrcv-Smirncv Test 
is a pair-wise test and The Krus kail- Wallis Test is a 
multi-sa I pie test tc determine if the underlying 
distriDutions are the same. 



1 . Kcuskall^Wallis Te st on Ra di al Errors 

A Kr uskall-Hallis Test was performed on the data to 
determine if the values were drawn from the same 
distribution. statistic (H) is generated ty merging all of 
the sample data, rank crdering the data and then grcupixig 
the data into their original sets and summing the ranks for 
each set. The Kruskall-Wallis statistic, H, is given ty 



H = (12/N(N + 1))(Rl/n + R2/m + R3/1 + Eh/p) - 3(N + 1) 
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Where : 





R1 = 


Sum 


cf 


ranks 


fcr 


Map 


A 








B2 = 


Sum 


cf 


ranks 


for 


Map 


E 








R3= 


Sura 


cf 


ranks 


fcr 


Map 


C 








E4= 


Sum 


cf 


ranks 


fcr 


Hap 


c 








n = 


tt = 


1 


= p 


96 












N = 


n ♦ 


12 


+ 1 + 


P = 


384 








The 


dis 


trib 


ut: 


ion of 


H is 


th 


en Chi-square 


with 


three degrees 


cf 


freed 


cm. 




















Ccrap 


arisons 


of 


the 


radial errors 


fcr 


the fcui maps 


for 


dayl 


ight 


a: 


nd night na 


ivigaticn trials 


cculd 


not reject 



the hypothesis that the radial errors for daylight cr the 
radial errors for night caoie from the same, underlying 
distribution . 



2. Kolmcgorov-S pi rnov Test of Sadisi 



a. Eaylight Navigation Trials 



Einperical Cumoiulati ve Distribution Functions 
were calculated for eacn map and plotted cn the same scale. 
This was accomplished by incrementing by five meters, 
counting the number of values which were egual to or less 
than that increment sum and dividing by the sample si 2 e of 
ninety-six. Using pair-wise calculations and an Alpha value 
of 0.10, the results were as follows. 

Let d=Max]Fl (X) -F2 (X) 1 
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Accept at Alpha=0.10 if d < 1. 22 (3grt (2/Sc) ) =0 . 176CS 



Jiap A 


vs Map E, d=0.06<0 . 17609 


I'iap A 


vs Map C, d=0.06<0 . 17609 


Map A 


vs Map D, d=0. 14<0 . 176GS 


Map B 


vs Map C, d=0.06<0 . 17609 


Map B 


vs Map D, d=0. 14<0. 17609 


Map C 


vs Map D, d=0. 12<0 . 176CS 



Ihe Kolincgc icv-Smirncv Test cannot reject the hypothesis 
that the underlying distributions are the saie. 
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VI. SIKULftTED ^STHIBUTTcns 



After ctserving the histograms of the actual radial 
errors and recognizing the fact that the radial errors viculd 
te positive values, the siaiulaticn of random variables which 
had values less than zero ware eliminated from 
consideratioc . The CDEC report used a Chi-Square Goodness 
of Fit Test and estimated the distribution of errors to be 
Exponential and was, therefore, included in the 
investigation. The Gamma distribution was investigated to 
try to explain the shape of the histograms and the larger 
number of values near zero. A lack of time prevented the 
investigation of possible combinations of distributions. 



A. SIMULAIICN OF THE EXPOHEKIIAL DIST5IBUIIC fJ 



1. General 



As Stated in Chapter II, CDEC performed Chi-Sguare 
Goodness of Eit tests to the radial errors ard stated that 
the exponential distribution provided the test fit. Since 
the Chi-Sguare Test is sensitive to how the date was 
sectioned and grouped at the tails of the distributions, 
this test is not very suitable. 
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2 . 



ISi S tatist ics lHves tiqat€ d 



Ey reading the radial errors into 
Library package, a histogram of rhe data was 
twenty-sii pertinent statistics were generated, 
coefficient of skewness, and coefficient cf 
three of the twenty-six statistics and were 
shewn belcfc. 



the HI5IF/g 
produced and 
The mean, 
kurtosis are 
computed as 



aean= X(i)/N 

Skewness= (third Central dement) /(Standard Eevation) 

Xurtosis= (Fourth Central Mcme nt) / (S tandard Devation) - 3 

The skewness and kurtosis of distributions are of interest 
because they have ranges .of values fer particular 
distributions. A range of values was obtained to compare 
with the real data by simulating the occurence of the random 
variable ever many replications. 



■3 • ijjb hod of S imu lati on 

The simulation procedure used was to calculate the estimate 
of -nhe mean for actual radial errors and to use this 
estimate as the parameter for an Exponential random variate 
generator. The random variate generator is described below 

X=Random Variate 

U = Uniform (0,1) random number 

M=Estimate of the Mean 

X=- (M) In (U) 
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Ninety-six of these Exponential (X) variates were generated 
tc simulate the sample size cf real navigation radial 
errors. These generated samples were then replicated one 
thousand times. 

4 . Ccmr ilation cf R epl ica tions 

Each sample of size 96 was read into a modified 
HISIf/Q sutroutine in crder to generate the sample 
statistics. The one thousand replications provided arrays 
for each statistic of size one thousand. The arrays fcr the 
mean, coefficient of skewness, coefficient of Kurtosis, 
Coefficient of Variation, Median, 50th Quantile, 75th 
Quantile and 90th Quantile were then put into the HISTF/G 
Liharary program. The output provided a sample distritution 
cf the statistics which were generated by an Exponential (M) 
distribution through which the actual data were compared. 
The tables cf these results are 'in Appendix E. 

E, SIMULAIICN OF THE GAMMA CISTEIBUTICN 



cf Parameters 

Ey inspecting the histogram generated by the 
experimental radial errors, the Gamma distribution appeared 
a logical possibility. To compute the scale and location 
parameters for the Gamma distribution, the calculated mean 
and variance were used as follows; 
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Wean=Sha pe* Scale 



Variance =Shape* (Scale) 

Ihis resulted in Alpfaa(Shape) and Beta (Scale) parameters 
calculated as follcws; 

Be ta= Variance/Mean 

aipha=Wean/B€ta 

Ihe calculated parameters of the distributicn are given in 
Table IV. 

TABLE IV - PAHAMETEFS OF GAMMA DISTEIEUTION 





DAYLIGHT TRIALS 


NIGHT TRIALS 


Map A; 


Alpha=0 .43901 


Alpha=0. 59214 




Bsta = 95 .99487 


Beta=78. 2295? 


Map E: 


Alpha=0 . 46086 


Alpha=0.8T735 




Seta = 96 .00516 


Beta=77. 17915 


Map C: 


Alpha=0. 59832 


Alpha=0. 53357 




Beta=78 . 94075 


Beta=102. 33755 


Map D : 


Alpha=0 .61320 


Alpba=1. 09553 




Beta =76. 92229 


Beta = 57. 5050 2 



2 • M et hod of Ga mja Simul aticn 

In order to generate sample sizes of ninety-six and 
replicate these sam-ples one thousand tines, a rejection 
method was utilized to generate Gamma variates. The Gamma 
random variate generator used was from the Directory of 
Routines in the International Mathematical and Statistical 
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Library (IM£L) at the X. E. Church Coaputer Center, Naval 
Postgraduate School, llcnterey, California. Since all IKSL 
routines are proprietary, the program cannot be included. 
Each of the samples of size ninety-six were then used as 
entries into the Modified KISTG/f Subroutine to generate 
sample sizes of one thousand for the twenty-six parameters 
for both the Daylight and Night Navigation errors. 



3 . CcE_g il^xo n of Re pli cati on s 



by the 
are in 



The compilation procedure was identical to that used 
Exponential simulation. Results of the compilation 
Appendix C. 



C. COMPAEISCN OF THE TWO SIMULATED DISTRIELIIONS TO THE 
ACTUAL DATA 



The comparisons below are made between the four map 
products for the parameters given. Replications cf tne 
parameters provided Normally distributed data for all tut 
the Kurtcsis and the Range. The support for a particular 
distribution was chosen by letting Alpha egual 0.1C and 
performing the standard Normal Test. Replications cf the 
simulated distributions did not provide Normally distributed 
Kurtosis or Range, sc standard application of statistical 
technigues was not appropriate. All Kurtcsis and Range 
parameters were checked to insure feasibility and not used 
to discriminate between the maps. Since both parameters for 
generating simulated Exponential the means of the 
distributions are not used as a discriminating parameter. 
The distribution sta tistics for the exponential and gamma 
distributions presented in tables V through XII are coded as 
fellows : 
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AV - The Actual Value. 

EMN - The Exgcnentiall j Generated Mean cf the Parameter. 
ESD - The Standard Deviation of the Exp Generated Mean. 
GaN - The Gamma Generated Mean of the Parameter. 

GSD - The Standard Deviation cf the Gamma Generated Mean. 

TABLE V - DISTRIBUTION STATISTICS - DAYLIGHT TRIALS - I'AE 



Parameter 


AV 


EMN 


ESD 


GMN 


GSE 


Mean 


h2. 14 


42.23 


4.233 


41.87 


6.523 


Std Dev 


63.60 


41 .91 


5.814 


62.29 


12.23 


Cceff Var 


1.509 


0. S93 


0.101 


1.488 


0.175 


Range 


318.7 


217.2 


53.. 66 


347. -7 


109. 2 


Skewness 


2.516 


1. 623 


0 .558 


2.632 


0.762 


Kurtosis 


4.257 


3.732 


8.7.27 


6.74 0 


6 .815 


50th Qnt 


19. 35 


29.55 


4.376 


17,24 


4 . 432 


75th Qnt 


38.97 


58.58 


7 .563 


54.23 


10.98 


90th Qnt 


115.8 


93.20 


12.03 


110.7 


21.44 



A 



The generated parameters for Map A shcfcn in Table V, 
support the Gamma distribution of the radial errors for this 
map. The Standard Deviation, Coefficient of Variation and 
all guantiles consistently reject the Expcner.tiai 
distribution and accept the Gamma distributicc. 
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TABLE VI 



CISTRiaUIICK STATISTICS 



DAYLIGHT TRIALS 



MAE 3 



Earamete r 


AV 


EKN 


ESD 


GMB 


GSD 


Mean 


44.25 


44. 17 


4.637 


41.87 


6.523 


Std Dev 


65. 18 


43.71 


6.441 


63.97 


12.24 


Coeff Var 


1.4 73 


0. 989 


0.098 


1. 448 


0. 176 


fiance 


356.8 


227.1 


56.11 


357.7 


116.5 


Skewness 


2.687 


1.S09 


0.572 


2.632 


0.762 


Kurtesis 


7.802 


4. 237 


3.949 


8.727 


6.740 


50th Qnt 


17. 73 


ao.96 


4 .420 


19.12 


4.795 


75th Qnt 


54. 15 


61.70 


7.94 1 


54.23 


10.98 


90th Qnt 


• 107.8 


97.74 


13.26 


110.7 


2 1.44 



The generated paraireters for Map 3, shcfcn in Tahle VI, 
again strongly support the Gamnia distribution of rhe radial 
errors for this map. The Standard Deviation, Coefficient of 
Variation and 50th Cuantile support a Gamma distribution. 
The rest of the parameters support either Gamma or 
Exponential distribution. The Skewness acxual value is 
within a normal range for both generated Exponential and 
Gamma. Eor the Experiential, the actual value is 1.535 
Standard Deviations frem the generated mean. For the Gamma, 
the actual value is 0.070 Standard Deviations frenr the 
generated mean. Clearly, the actual value is more 
"ccDif ortable" in the Gamma distribution, but for any 
reasonable Alpha value (<0.10) the hypothesis that the 
actual value is from the Exponential distribution canrot be 
rejected. The 75th and 90th Quantiles egually support 
either a Gamma or Exponential distribution. 
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lAELE VII - EISTEIBUTICN STAIISTICS - DAYIIGHI TRIALS - MAf 
C 



Parameter 


AV 


EKN 




Heati 


47.23 


47.40 


4 


Std Dev 


61.06 


46.87 


6 


Cceff Var 


1.293 


0. 988 


0 


fiange 


273.6 


241.7 


5 


Skewness 


2.062 


1.782 


0 


Kurtosis 


3.654 


3. 967 


3 


5Cth Qnt 


2 2.34 


53.24 


4 


75th Qnt 


51.92 


66.03 


8 


90th Qnt 


138.7 


105.1 


1 



The generated faraireters 
support the Gamma distributi 
map. The Standard Deviation 
Quantile and 90th Quantile 
The Skewness and 75th Quanti 
Gamma distribution. For t 
0.555 Standard Deviations f 
0.233 Standard Deviations fr 
Quantile, the actual value i 
the Exponential mean and 1. 
Gamma mean. 



ESD GliN 


GSE 




.936 47.45 


6.146 




.666 59.93 


10.26 




.095 1.262 


0. 129 




7.54 318.2 


108. 4 




.506 2.216 


0.618 




.276 6,106 


4.74? 




.902 25.42 


5.193 




.254 64.76 


10. 18 




4.09 118.4. 


19.54 




for Map C, 


shown in Table 


VII, 


on of the radial errors for 


this 


, Coefficient of Variation, 


5Cth 


support the 


Gamma distribution. 


le support e 


ither Exponential of 


he Skewness, 


the actual valu 


e is 


rom the Exponential mean 


and 


cm the Gamma 


mean. For the 


75th 


s 1.709 Standard Deviations 


from 


261 Standard 


Deviations from 


t he 
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lABlE VIII - CISTRIBUTION STMISTICS - DAYLIGfiT THIAIS - MAP 
C 



Parameter 


AV 


EKN 


ESD 


GMN 


GS£ 


Mean 


47.17 


47.22 


4.651 


47 .23 


6.381 


Std Dev 


60.24 


46.51 


6.53 7 


59.35 


1C. 18 


Cceff Var 


1.277 


0.985 


0.095 


1.25 6 


0 .131 


Range 


281.0 


239.9 


59.55 


318.4 


88.99 


Skewness 


1.831 


1.787 


0.534 


2.256 


0.656 


Kurtosis 


2.963 


4.C50 


3.482 


6.438 


5.027 


50th Qnt 


17. 46 


33.12 


4 .751 


25.67 


5.155 


75th Qnt 


72.01 


65 .74 


8 .081 


64.54 


10.66 


90th Qnt 


126.6 


104.4 


12.99 


116.7 


19.60 



The generated paraceters for Map D, shckc in Table VIII, 
support the Gamma distribution of the radial errors for this 
map. The Standard Deviation, Coefficient of Variation and 
50th Quantile support the Gamma distribution. The remaining 
parameters support either a Gamma or Exponential 
distribution. Except for the Skewness, the Gamma 
distribution fit is better than the Exponential 
distribution. The Exponential simulation provides a tetter 
fit of the actual Skewness. 
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lABlE IX - CISTfilBUTIOK STATISTICS - NIGHT TBIALS - MAE A 



Earamete r 


AV 


EMN 


ESD 


GMN 


GS£ 


Mean 


46.32 


46. 42 


4.653 


46.33 


5 .944 


Std Dev 


60.20 


.46.06 


6 .391 


59.65 


1.018 


Cceff Var 


1.300 


0. 993 


0.100 


1.287 


0.141 


Range 


298.0 


238.7 


59,81 


32^. ^ 


92. 13 


Skewness 


2.268 


1.632 


0.558 


2.322 


0.669 


Kurtosis 


5.283 


4.257 


3.732 


6.860 


5.53C 


50th Qnt 


19. 26 


32.48 


4.81C 


24.50 


4 .929 


75th Qnt 


61.48 


64.39 


8.315 


62.78 


10.37 


SOth Qnt 


109.0 


102.4 


1 3.23 


115.6 


18.63 



The generated parameters for Map A, shc«in in Table IX, 
support the Gamma distribution cf the radial errors for this 
map. The Standard Deviation, Coefficient cf variation and 
50th Quantile support the Gamma distribution. All of the 
remaining generated parameters support either the 
Exponential or Gamma distribution. 
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TABLE X - DISIBIBUTION STATISTICS - MIGHT lEIALS - ilAP E 



Parameter 


AV 


EKN 


ESD 


GMN 


GSE 


Mean 


63.08 


62 .98 


6.612 


62.97 


7.138 


Std Dev 


69.78 


62.33 


9.184 


68.62 


10.47 


Coeff Var 


1.106 


0.969 


0.098 


1.088 


0.114 


Bange 


463. 1 


323.8 


82.85 


360.4 


96.13 


Skewness 


2.773 


1.809 


0.572 


2.322 


0.689 


Kurtosis 


11. 33 


4. 237 


3.949 


5.074 


4.629 


50th Qnt 


45.79 


44.15 


6.302 


40.62 


6.758 


75th Qnt 


87.23 


87.98 


11.32 


37.77 


12.12 


SOth Qnt 


134. -4 


139.4 


18.93 


145.9 


21 .03 



The generated paraneters for Map B, sbciir in Table X, 
support either the Gaaica or Exponential distribution of the 
radial errors for this map. No statistical discriminator 
can be found to determine which distrituticn is correct. 
This fact is not surprising, since the snape parameter for 
the Gamma distribution was C.817, close to the shape 
parameter for the Gaamia to become an Exponential (Shape = 
1 . 000 ) . 
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lAflLf XI 



DISIHIBUTION STATISTICS 



NIGHT T5IALS 



MAE C 



Earameter 


AV 


EKN 


ESD 


GMN 


GSI 


Mean 


54.60 


54.79 


5 .708 


55.15 


7.427 


Std Dev 


74.75 


54.18 


7.706 


73.9 0 


13.04 


Coeff Var 


1.3 69 


0. S89 


0 .095 


1.340 


0. 151 


Bange 


557.4 


279.4 


66.41 


397.8 


117.6 


Skewness 


3.818 


1. 782 


0.572 


2. 352 


0.6SS 


Kurtosis 


3.654 


3. S76 


3.27 6 


6.106 


4.747 


50th Qnt 


33.58 


38.42 


5 .667 


27.01 


5.95C 


■ 75 th Qnt 


62.60 


76 .33 


9.542 


74.20 


12.76 


SOth Qnt 


135. 1 


121 .5 


16.29 


141 . 0 


24.61 



The gen€rated parameters for Map C, shckn in Table XI, 
support the Gamma distribution of the radial errors for this 
map. The Standard Deviation, Coefficient of Variation and 
Skewness support the Gamma distribution. The Quantiles 
support either an Exponential or Gamma distribution. 
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lABLE XII 



CISTRIBUTICN STATISTICS 



NIGHT TRIALS 



MAP D 



Parameter 


AV 


EMN 


ESD 


GMN 


GSC 


Mean 


63.00 


63.07 


6.212 


40.20 


3. 962 


Std Dev 


60.19 


62.12 


8.731 


37.92 


5.153 


Coeff Var 


0.985 


0. S85 


0.095 


0.944 


0.090 


Range 


302. 3 


320.5 


79.53 


194.2 


44.55 


Skewness 


1.490 


1.787 


0.534 


1.698 


0.507 


Kurtosis 


2. 253 


4. 050 


3.482 


3. 652 


3.214 


5Gth Qnt 


50. 26 


44 .23 


6 .346 


29.14 


4.069 


75th Qnt 


86.79 


87.81 


10.79 


56.11 


6.752 


90th Qnt 


161.0 


139. 5 


17.35 


87.18 


11.22 



The generated parameters for Map B, shc^n in Tatle XII, 
support the Exponential distribution of radial errors for 
this map. The Standard Deviation and all guantiles support 
the Exponential distribution. Only the Coefficient of 
Variation and Skewness support either the Gamma or 
Exponential distribution. The resulting support for the 
Exponential distriburion is not unexpected given the mean 
and standard deviation relationship in the actual values. 
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VII 



£'CNCIii^SIONS 



A. IHE CISIEIBOTICN OF RADIAL ERRORS 



The siiiulated gatuna distributicn provid 
of the parameters tc the radial errors. 
Exponential distribution is a special ca 
distribution/ it was net surprising that t 
distribution provided good fits to the rad 
CCEC performed the Chi-Sguare Goodness of Fit 
night navigation trials, is probably best 
the Expccential distribution. For all ctfc 
maps, the errors are best approximated 
distribution. 



ed the best fit 
Since the 
se of the Gamma 
he Exponential 
ial errors when 
Test. Map D, 
approximated by 
er trials and 
by the Gamma 



E. 



COMRAEISCN OF THE 



Maps BY QUANTILES 



Sine 
between 
was mad 
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"* • ' !th Qua ntile 

Map A's real quantiles were statistically tetter 
than the ether three maps during the daylight trials and 
tetter than all but Map D during the night trials. “Thus, 
when flap & was used, 75 percent of the values were 
significantly better than the other maps. 



2- The SCth Qua nt ile 

No statistical difference exists between the maps at 
the 90th guartile during the daylight trials. The night 
trials produced a difference only in that Map D was the 
statistically worst map. Based upon this statistic, Mac D 
was eliminated as the worst alternative airccg the four map 
products. 

C. COflPABISCN OF THE MAPS BY AREAS CF THE EE50RS 

An alternative method of ranging map products involves 
returning to the X and Y values. Because of the methcc used 
to generate the X and Y values, a coinmcn checkpoint was 
designated and the navigation errors can he placed on a 
single axis for each map product and trial. Observation of 
these errors on a common plot leads to another method of 
ranking. Since the goal of MAPPEO III was to determine the 
performance cf the map products, one measure could he tc 
rank the maps by the size cf the areas generated by a 
percentage of the results. Ssmemberino that positive and 
negative values were available for X and Y values, the lower 
percentiles would , in general, reflect negative values and 
the higher percentiles would reflect positive values. By 
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calculating the absolute distance between lower and higher 
percentiles, a percentage of values will fall intc that 
absolute distance. By doing this for both X and Y values 
and then multiplying the absolute distarce for X by the 
absolute values for Y, an area which includes an 
approximation to a percentage of the (X,Y) values in that 
area results. The size of the area then beccmes an analysis 
tool to rank the maps. An analogy exists with the 
percentiles given in previous sections for radial errors. 
Table XIII provides the results of calculating the areas of 
50-percent (75th Percentile minus the 25th Eercentile) and 
80-percent (90th Percentile minus the 10th Percentile) 
areas. The figures are in meters squared and the rank is 
given with tne smallest area being ranked 1 and the largest 
area being ranked 4. Fortunately, differences did net occur 
in the rankings f cr 5C-percent and 80-percent. Map A ranks 
first for both day and night trials as it did using radial 
errors. 
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TfiBLE XIII - EANKISGS OF HAPS EY A3EAS 
DAYXIGHT TRIALS 



HAP 


50PCT 


80PCI 


RAN 


A 


304 


1754 


1 


B 


34 2 


6293 . 


2 


C 


407 


7735 


4 


D 


407 


7735 


3 




NIGHT 


TRIALS 




A 


708 


8284 


1 


B 


2365 


13640 


3 


C 


1300 


9477 


2 


D .. 


2336 


15300 


4 



The A NOVA technigue used by CDEC was af f ropriate f cr the 
distribu ticQ of the error data. However, a ranking cf the 
maps during navigational trials was not possible. Through 
the other techniques employed , a ranking is possible for the 
navigation trials. Hap A, the map currently in use, is the 
best of the alternatives presented. 

C. APPLICATION OF THE RESULTS 

The sample size for the MAPPfO III experiment was not 
large enough to reach a definite conclusion on the data. 
With the date available, the conclusion that the radial 
errors experienced during navigational trials is Gamma 
distributed is valid. This knowledge can be used during 
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future (and past) field experiments with wap products to 
perform parametric tests on the estimated paraneters 
generated by the actual data. The parameters for the Gamma 
distribution should be calculated as shown tc provide the 
test results. Use cf the Fortran routines supplied will 
permit the simulation of the distribution to check the 
actual data. Once this simulation is performed and the 
parameters meet the desired criteria, parametric tests can 
te performed to evaluate the data. 
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APPENDIX A 



KHUSKAL-WALLIE TEST BATA 

NOTE: The following tables provide the test data froc the 
Kruskal- Siiallis Test. The percents given are produced from 
the Cumulative Distribution Function which has values from 0 
to 1.00. 
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DAYllGHT TRIALS 



(0,5) 


:iAP 

A 

ECT 

20 


(0,10) 


35 


(0,15) 


42 


(0,20) 


53 


(0,25) 


61 


(0,30) 


67 


(0,35) 


69 


(0,40) 


72 


(0,45) 


75 


o 

o 


76 


(0,55) 


76 


(0,50) 


77 


(0,65) 


78 


(0,70) 


79 


(0,75) 


82 


(0,80) 


86 


(0,85) 


89 


(0,90) 


90 


(0,95) 


90 


(0,100) 


92 



HAP 


MAP 




C 


D 


MAXIMUM 


PCT 


PCT 


DIFFEHHN 


14 


10 


10 


34 


22 


14 


48 


40 


08 


56 


49 


07 


64 


53 


1 1 


66 


58 


09 


67 


61 


08 


68 


63 


09 


71 


66 


09 


74 


69 


06 


78 


71 


09 


79 


73 


07 


79 


74 


08 


80 


75 


08 


80 


76 


07 


82 


77 


09 


83 


77 


12 


85 


73 


12 


86 


78 


12 


88 


78 


14 



«AP 

B 

PCT 

15 

36 

48 

54 

61 

66 

68 

71 

74 

77 

80 

80 

82 

83 

83 

83 

'85 

85 

86 

86 



50 





BAP 

A 

Ecr 


(0,105) 


92 


(0,110) 


92 


(0,115) 


92 


(0,120) 


94 


(0, 125) 


94 


(0,130) 


94 


(0,135) 


94 


(0,140) 


94 


(0,145) 


95 


(0,150) 


95 


(0,155) 


95 


(0,160) 


95 


(0, 165) 


96 


(0,170) 


96 


(0, 175) 


96 


(0,180) 


96 


(0,185) 


96 


(0,190) 


96 


(0,195) 


96 


(0,200) 


96 


(0,>200) 


100 



BA P 


HAP 




C 


D 


MAXIHUd 


PCI 


FCT 


DIfFEEHN 


89 


78 


14 


90 


81 


11 


90 


83 


09 


90 


84 


10 


90 


85 


09 


90 


85 


09 


90 


85 


09 


90 


88 


06 


90 


90 


05 


90 


90 


05 


90 


91 


05 


90 


93 


05 


91 


93 


05 


91 


93 


05 


92 


S3 


04 


92 


93 


04 


92 


94 


04 


92 


94 


06 


93 


94 


05 


93 


94 


05 


100 


100 


00 



HAF 

B 

FCT 

86 

86 

86 

89 

8S 

91 

92 

93 

93 

93 

94 

94 

94 

94 

94 

95 

96 

98 

98 

98 

100 



51 



NIGHI TRIALS 





MAP 

A 

ECT 


(0,5) 


14 


o 

o 


28 


(0,15) 


34 


(0,20) 


39 


(0,25) 


46 


(0,30) 


48 


(0,35) 


50 


(0,40) 


57 


(0,45) 


60 


(0,50) 


60 


(0,55) 


64 


(0,60) 


67 


(0,65) 


70 


(0,70) 


75 


(0,75) 


77 


(0,80) 


76 


(0,85) 


79 


(C,90) 


82 


(0/95) 


82 


(0, 100) 


84 



MAP 


MAP 




C 


D 


MAXIMUM 


PCT 


PCT 


DIPPERLN 


10 


09 


06 


22 


24 


08 


30 


35 


05 


41 


39 


09 


43 


41 


10 


45 


47 


08 


48 


50 


05 


48 


55 


09 


. 49 


58 


11 


55 


64 


09 


65 


69 


11 


69 


73 


13 


70 


75 


09 


73 


77 


10 


75 


77 


06 


76 


80 


05 


77 


80 


04 


80 


60 


04 


83 


81 


04 


85 


81 


04 



J«AP 

B 

ECT 

08 

20 

31 

32 

36 

40 

45 

4-8 

•5 4 

56 

58 

63 

66 

67 

71 

75 

76 

78 

79 

85 
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MAP 

A 

PCT 



(0,105) 


86 


(0,110) 


86 


• (0,115) 


86 


(0,120) 


86 


(0,125) 


88 


(0,130) 


88 


(0,135) 


88 


(0,140) 


88 


(0,145) 


89 


(0, 150) 


89 


(0,155) 


90 


(0,160) 


91 


(0,165) 


91 


(0,170) 


91 


(0,175) 


91 


(0, 180) 


91 


(0,185) 


92 


(C, 190) 


93 


(0,195) 


95 


(0,200) 


96 


(0,>200) 


100 



MAP 


MAP 




C 


D 


MAXIMUM 


PCI 


PCI 


DIFFEBEM 


89 


82 


07 


89 


82 


07 


89 


83 


07 


89 


85 


05 


90 


35 


06 


91 


85 


06 


91 


85 


06 


91 


86 


07 


92 


88 


07 


92 


88 


07 


94 


88 


08 


94 


88 


08 


96 


89 


09 


96 


89 


09 


96 


89 


09 


96 


89 


10 


96 


90 


09 


97 


92 


08 


97 


93 


C7 


97 


93 


07 


100 


100 


00 



MAP 

3 

PCI 

89 

89 

90 

90 

91 

91 

9 1 

93 

95 

95 

96 

96 

98 

98 

98 

99 

99 

100 

100 

100 

100 
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APPENDIX B 



PABANEIEfiS GEiiEBAIEC FOE THE EXPONZNTIAI DISTRIEUIICN 

NOTE: The results given are rrcm one thousand replications 
of 96 Expcaential generated random variates. All values 
except fcr skewness, kurtosis and coefficient of variation 
are in meters. The skewness, kurtosis and coefficient of 
variation are diaensiorJ.ass quantities. 
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RESULTS EOB DAYLIGHT NAVIGATION 



A. 



1 . 


Generated 


MeanSi 1000 Replications 








MAP A 


MAP B 


MAE C 


MAP 


Mean 




42.23 


44.10 


47.40 


47. 


Standard 


Deviation 


4, 23 3 


4.637 


4 .938 


4.6 


Skewness 




0.01 7 


0.294 


0.239 


0.0 


Kurtosis 




0. 13 1 


0.526 


0.367 


-0. 


tiinimum 




29.13 


31.00 


31.02 


3 ^ . 


Maximum 


• 


57.12 


65.47 


71.22 


59. 


2. 


Generated 


Standard 


Deviation, 


J.0d0 Replica 


ticn 






MAP A 


MAP B 


MAP C 


MAP 


Mean 




4 1.9 1 


43.71 


46.67 


46. 


Standard 


• Deviation 


5.81 4 


6.441 


6.666 


6.5 


Skewness 




0.437 


0.504 


0.563 


0.4 


Kurtosis 




0.36 9 


0.678 


0. 8SC 


0. 1 


Minimum 




24.68 


27.32 


28.40 


30. 


Maximum 




64.03 


73.30 


75. SS 


71, 



D 

22 

51 

35 

16 

i)0 

27 

c 

D 

51 

37 

20 

20 

C3 

S2 
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i, 
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• Generated Skew ness . lOCO Replica t i cns 







MAP A 


HAP B 


MAP C 


MAP D 


Mean 




1 .82 3 


1 .809 


1.782 


1.787 


standard 


Deviation 


0.55 8 


0 .572 


0.5C6 


0.534 


Skewness 




1 .08 8 


1.316 


1.C63 


1. CG4 


Kurtosis 




1 .95 8 


2.783 


2.270 


1.682 


Kinimu.m 




0.72 7 


0.632 


0.728 


0.717 


M a ^ 1 [D u Qi 




4.75 2 


4.822 


4.354 


4.654 


4. G 


enerated 


Kurtosis , 


1 00 ( Repli 


cations 








MAP A 


MAP B 


MAP C 


MAP D 


Mean 




4 .25 7 


4.237 


3.976 


4.050 


Standard 


Deviation 


3.732 


3.949 


3.278 


3.482 


Skewness 




2 .22 0 


2.485 


2.259 


1.983 


Kurtosis 




7.77 1 


9.171 


8.565 


5.742 


MiniQum 




-0.42 


-0.81 


-C.42 


-0.60 


Maxiuuiii 




32.27 


33.03 


27. 13 


26.82 
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5- Gaijer^ed Coeffici en t cf Var iati on^ 



Eeolicat 


icas 














MAP A 


MAP B 


MAP C 


MAP D 


Mean 




0.99 3 


0.989 


0. see 


0.985 


Standard 


De V iation 


0.101 


0.098 


0.095 


0.095 


Skewness 




0.629 


0.467 


0.481 


0.546 


Kurtosis 




0.721 


0.532 


0.393 


1.132 


Mi aim urn 




0.75 2 


C.689 


0.714 


0.731 


Ma ximum 




1 ,44 4 


1 .410 


1 . 345 


1.486 


6. 


Generated 


5Cth Ouantiie 


!x 


Beplications 








MAP A 


MAP B 


MAE C 


MAP D 


Mean 




29.55 


30,96 


33.24 


33. 12 


Standard 


De V iation 


4 .37 6 


4.420 


4.902 


4.751 


Sk-ewness 




0.299 


C.257 


0.2S1 


0.3 13 


Kurtosis 




0.102 


0,374 


0. C84 


0.172 


Minimum 




1£.18 


18.07 


20.20 


20. 55 


Maximum 




44.07 


47.76 


54.09 


50.29 
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• <3€ner at e4 75th Q uan t ile, 1000 Replicat ions 







MAP A 


HAP B 


MAP C 


MAP D 


Mean 




58.58 


61.70 


66 . C3 


65.74 


Standard 


Deviation 


7.56 3 


7.941 


8.254 


8.081 


Skewness 




0.287 


C.312 


0.240 


0. 1S6 


Kurtosis 




0.30 8 


0.312 


0.2C7 


0.064 


Minimum 




37.21 


41.37 


4 0.31 


42. 80 


Maximum 




88.21 


93.81 


102.0 


94.57 


8. 


Generated 


SGth Quantile 


\jL 1QQQ 


Replications 








MAP A 


MAP. 3 


MAP C 


MAP E 


Mean 




93.20 


97.74 


ICS. 1 


104.4 


Standard 


Deviation 


12.03 


13.28 


14. OS 


12. S9 


Skewness 




0.291 


0.369 


0.352 


0. 25S 


Kurtosis 




-0.04 


-0.U9 


0.616 


C.2S0 


Minimum 




61.97 


67.00 


6 8.06 


62. 43 


Maximum 




133. 7 


142.2 


181.2 


156.0 
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9. Qsneratsd 


pan^SjE^ 1000 


Replicati 


pns 




MA? A 


MAP B 


MAE C 


Mean 


217. 2 


227.1 


241.7 


Standard Deviation 


53.66 


58. 11 


57.54 


Skewness 


1.002 


1.163 


1.066 


Kurtosis 


1 .33 4 


2.185 


2.215 


Minimum 


96.69 


1 17.5 


120.7 


Ma:cimum 


47C. 0 


518.2 


527.7 



MAE B 
239 .9 
59.55 
0.9'?4 
1.462 
128.5 
544.2 
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c. RESULTS FOR NIGHT NAVIGATION 



”• • ed HeanSi lOOC Rep licat ions 





HAP A 


MAP B 


MAP C 


MAP D 


Mean 


46.42 


62.98 


54.79 


63. C7 


Standard Deviation 


4 .65 3 


6 .612 


5 .7C8 


6.2 12 


Skewness 


O '. 07 1 


0.294 


0.239 


0.035 


Kurtosis 


0.131 


0.526 


0.367 


-0. 16 


Minimum 


32.02 


44.21 


35.66 


43.27 


Uaxlmum 


62.77 


93.35 


82.33 


79.15 


2.- Generated 


Stan dard 


Deviation, 


1000 Replica 


ticne 





MAP A 


MAP B 


MAP C 


MAP D 


Mean 


4 6.06 


62.33 


54. 16 


62. 12 


Standard Deviation 


6 .391 


9.184 


7. 706 


8.731 


Skewness 


6-437 


0.504 


0.563 


C. 420 


Kurtosis 


0.36 9 


0,678 


0. 8S0 


0- 120 


Miniaiuai 


27,13 


38.96 


32.83 


4-0.11 


MaxiiDuni 


7G.38 


1C4.5 


87.85 


96. C5 
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3. 


^§Sjrjted Skgwness « 


1552 Beplica~!;icc£ 






MiiP A 


MAP B 


MAP C 


MAP D 


il€an 


1 .82 3 


1.809 


1 .782 


1.787 


Standard 


Deviation 0.558 


0.572 


0.5C6 


0.534 


Skewness 


1 .08 8 


1.316 


1 . 063 


1 .OCC 


Kurtosis 


1 .95 8 


2 .783 


2.270 


1.682 


Kinioiuai 


0.72 7 


0.632 


0.728 


0.717 


I'l a X i !D u m 


4.75 2 


4.822 


4.354 


4.654 


4. 


Generated Kurtosis* 


1000 Replica 


t ions 






MAP A 


MAP E 


HAP C 


MAP D 


Mean 


4.25 7 


4 .237 


3.S76 


4.050 


Standard 


Deviation 3.732 


3 .949 


3 .278 


3.482 


Skewness 


2.22 0 


2.485 


2.259 


1.963 


Kurtosis 


7.771 


9.171 


8.565 


5.742 


Minimum 


-0.42 


-0.81 


-0.42 


-0.59 


Maximum 


32.27 


33. 03 


27. 13 


26.62 
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Variation 



1000 



5* gen era ted Coeffici^t of 
Seplicatiens 





MAP A 


MAP B 


MAP C 


MAP D 


Mean 


0.99 3 


0.989 


0.569 


0.965 


Standard Deviation 


0.100 


0.098 


0.C95 


0.055 


Skewness 


0.62 9 


0.467 


0.481 


0.546 


Kurtosis 


0.721 


0.532 


0 .353 


1.132 


Minimum 


0.75 2 


C.689 


0.714 


0.731 


Maximum 


1 .444 


1 .410 


1.345 


1.486 


Generated 


5^th Quantile 


•X 


Replications 




• 


. MAP A 


MAP B 


MAP C 


HAP D 


Mean 


32. 4S 


44. 15 


36.42 


44-23 


standard Deviation 


4 .81 0 


6.302 


5.667 


6.346 


Skewness 


0.299 


C.310 


0. 2S1 


0.3 13 


Kurtosis 


0.102 


0.115 


0 . C64 


0.172 


Minimum 


15.98 


20.05 


23.35 


27.45 


Maximum 


48.44 


49.61 


62.53 


67. 17 
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7. Ggneratad 


75th Ouantil 


1000 


Replicati 


ons 






MAP A 


MAP B 


MAE C 




HAP D 


Mean 


64.39 


87.98 


76.33 




87.61 


Standard Deviation 


8.313 


11.32 


9.542 




10.79 


Skewness 


0.28 7 


0.312 


0.240 




0.156 


Kurtosis 


0.308 


0.312 


0. 2C7 




0.064 


Minimuni 


4C.90 


58.99 


46. 6C 




57. 17 


Maximuai 


96.96 


133.8 


117.9 




126.3 


Generated 


J^^ii Quant il 


§x JOOO 


Reglicati 


ons 






HAP A 


MAP B 


MAE C 




MAE C 


Mean 


102. 4 


139.4 


121.5 




139.5 


Standard Deviation 


13.23 


18.93 


16.29 




17.35 


Skewness 


0.291 


0.369 


0.3S2 




0.255 


Kurtosis 


-C.04 


-0.09 


0.616 




0.250 


Miniamm 


66.12 


95.54 


78.67 




83.38 


Maximum 


146. 9 


202.7 


209.5 




208.3 
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9 • Geae r a t€ d Hangs, 1000 Re plicat ions 



MAP A 

Mean 238. 7 

Standard Deviation 5fc.98 
Skewness 1.002 

Kurtosis 1.334 

Minitnum 1C8.5 

Maximum 5 16. 6 



MAP B 


MAP C 


MAE C 


3 23 .' 8 


379.4 


320.5 


82.85 


6 6.41 


79.53 


1 .168 


1.066 


0.974 


2-185 


2.215 


1.462 


167.6 


139.6 


171.6 


738.9 


6 10.0 


726.8 
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APPENDIX C 



PAEAJIETDRS GENERATED FOR THE GAaMA DISTRIBUTION 

NOTE: The results giver are from one thousand replications 
cf 96 Gamma generated random variates. All values except 
for skewness, kurtcsis and coefficient of variation are in 
meters. The skewness, kurtosis and coefficient of variation 
are dime nsicnless guantities. 
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RESULTS EOfi DAYLIGHT NAVIGATION 



1 . 


Generated Means, 1000 


Replications 






MAP A 


MAP B 


HAP C 


MAE D 


Mean 


4 1.67 


44.16 


47.45 


47.23 


Standard 


Ceviation 6.523 


6 .470 


6.146 


6.381 


Skewness 


0 . 25 1 


C.460 


0. 183 


0.350 


Kurtosis 


0.03 5 


0.382 


0.116 


0.253 


Miniaium 


25.13 


26.65 


28.52 


28.84 


Maximuni 


6 1 .36 


7C. 18 


69.43 


69. 64 


2. 


Generated Standard Dev 


ia tion , 


JOOO I'gplications 




MAP A 


MAP 3 


MAP C 


MAP D 


Mean 


62.29 


63.97 


59. S3 


59.55 


Standard 


Ceviation 12.23 


12.24 


10.26 


1C. 18 


Skewness 


0.45 2 


0 .657 


0 .5S6 


0.5C4 


Kurtosis 


0. 17 6 


C.844 


0.7C3 


C. 556 


M 1 n i 01 u IS 


27.92 


33.46 


35.33 


34.72 


Ma xiaium 


107. 2 


117.2 


107.6 


104.7 
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3. Gjjieiated Sk ew ness , 100 0 Repli c a t iccs 





a^P A 


MAP B 


MAE C 


MAE D 


Mean 


2.63 2 


2.579 


2-216 


2.256 


Standara Eeviation 


0.78 2 


0.806 


0.618 


0.656 


Skewness 


1.196 


1.374 


1 . 1 16 


1.0 17 


Kurtosis 


2.009 


2.398 


1 .899 


1.573 


Minimua 


1 .170 


0.983 


0.876 


1.012 


Maxiflium 


6.426 


6 .344 


4 .853 


5.6C8 


*♦ • Gene rj ted 


Kurtcsis . 


1000 Reel 


ications 






-MAP A • 


MAP B 


MAE C 


MAE D 


Mean 


8.727 


8.538 


6. 1C6 


6.458 


Standard Deviation 


6 .74 0 


6.998 


4.747 


5.027 


Skewness 


2. 112 


2.156 


2.115 


1.920 


Kurtosis 


5.994 


5.821 


6. 185 


5.398 


MinimuiB 


0 .37 9 


-0.26 


-0.42 


-0. C5 


Maximum 


50.22 


4 9.38 


3 3.25 


39.62 
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5 



of 



Vari atio n 



100 ^ 



Gare rate d Coeff icien t 
E€£l leaf lens 





MAP A 


MAP B 


MAP C 


MAP B 


Mean 


1 .48 8 


1.448 


1.262 


1.256 


Standard Deviation 


0.175 


0.178 


0 . 129 


0.131 


Skewness 


0.59 4 


0.833 


0.4C7 


0.417 


Kurtosis 


0.590 


1 .060 


0.087 


0.6C4 


Mininum 


1 .06 4 


1.072 


0.53S 


0.861 


Maximuni 


2. 13 9 


2.165 


1.754 


1.926 



6. 


Generated 5 Ctn Quant 


il6i 1000 


Replications 






M AP ■ A 


MAP B 


MAP C 


MAP D 


Mean 


17.24 


19.12 


25.42 


25. 67 


Standard 


Deviation 4.432 


4.795 


5.193 


5. 155 


Skewness 


0.718 


0.576 


0.365 


0.466 


Kurtosis 


0.88 0 


0.467 


O.OC1 


0.555 


Minimum 


7.856 


7 .835 


12. 6C 


11.78 


Maximum 


39.58 


38. 74 


45. 4C 


46. S7 
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7. Generated 


75th Quantile 


^ 1000 


Replications 






MAP A 


MAP B 


MAP C 


MAP D 


Mean 


54.23 


58.09 


64.76 


64.54 


Standard Deviation 


10.98 


10.99 


10.13 


10. 66 


Skewness 


0 .56 7 


0.500 


0 . 187 


0.4 63 


Kurtosis 


0.499 


0.337 


-0.12 


0.466 


Minimum 


• 29.66 


29.46 


34.73 


37.29 


Maximum 


100. 6 


98.50 


94.41 


109.5 


Generated 


90th Quantile 


, 1000 


Replications 







MAP A . 


HAP B 


MAP C 


MAP D 


Mean 


1 10. 7 


1 15.5 


118.4 


116.7 


Standard Deviation 


21.44 


20.86 


19.54 


19.60 


Skewness 


0.35 6 


0.555 


0 .484 


0.439 


Kurtosis 


0 .56 1 


0.457 


0.598 


0.574 


Minimum 


54.49 


67.09 


66.96 


66 . 83 


Maximum 


184. 6 


208.7 


200.1 


192 .8 
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9. 


generated 


T QO C 


Replication 


3 








MAf A 


MAP B 


MAP C 


MAP D 


Mean 




347. 7 


357.7 


318.2 


318 .4 


Stacdara 


Ce viation 


109. 2 


1 16.5 


87.98 


88.99 


Skewness 




1. 16 9 


1 .306 


1 . lie 


1.021 


Kurtosis 




2.233 


2.682 


1.831 


1.797 


ISinimum 




122. 4 


139.7 


153.2 


148.7 


riaxiaiuin 




9Ce. 6 


948.3 


732.1 


821 .4 
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RESULTS FOB NIGHT IIAVIGATICN 



B . 



Ms ans JOOO R e Flicati c ns 





HAP A 


MAP B 


MAP C 


MAP E 


Mean 


46.33 


62.97 


55.15 


40.20 


Standard 


Eeviaticn 5.S44 


7 .138 


7.427 


3.982 


Skewness 


0 . 16 4 


0.201 


0.219 


0.159 


Kurtosis 


0.02 9 


-0.02 


0. 188 


-C. C2 


Miniauat 


2t. 18 


42.92 


34.40 


28. 46 


Ha xiaua 


68.54 


89.24 


84.52 


52.91 


2. G 


eneiated Standard 


Deviation , 


iMO Beclications 




HAP A 


MAP B 


HAP C 


MAP C 


Kean 


59.65 


68.62 


73. SC 


37.92 


Standard 


Eeviation 1C. 18 


1 0.47 


13.04 


5.153 


Skewness 


0.555 


0 .499 


0.64 1 


0.420 


Kurtosis 


0.66 2 


0.307 


0.811 


0.229 


Kininum 


36.17 


42.32 


4C.87 


25. 64 


Kaxiinuffl 


1 11. 8 


1C6.3 


127.4 


57.79 
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3. 


Generated 


Skewness , 


JOOO Replies 


tiens 








HAP A 


MAP B 


MAE C 


MAP B 


Mean 




2.32 2 


1.979 


2.352 


1.698 


Standard 


Ee viation 


0 .689 


0.627 


0.699 


0.507 


Skewness 




1 .23 5 


1 .426 


1.275 


1.004 


Kurtosis 




2.024 


3.485 


2 .706 


1.772 


Minimuni 




1.071 


0.867 


0.946 


0.6C6 


Maxim um 




5.46 5 


5.201 


6.611 


4, 127 


4. 


Generated 


Kurtosis , 


JOOO Replications 








MAP A 


MAP B 


MAP C 


MAP D 


Mean 




6 . 85 0 


5.074 


6.932 


3.652 


Standard 


Deviation 


5.53 0 


4.629 


5.570 


3. 214 


Skewness 




2.11 5 


2.670 


2.273 


2.136 


Kurtosis 




5.707 


10.34 


8. 104 


7.146 


Minimum 




0.097 


-0.20 


-0.13 


-0.63 


Maximum 




37.64 


35.24 


52.32 


25.65 
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i 











5 • Generated Coeffi cien t cf Var iati on ^ 

SS£iicat icns 





MAP A 


MAP B 


MAP C 


MAP D 


Mean 


1 .23 7 


1 .088 


1 .340 


0.944 


Standard Deviation 


0. 141 


0.114 


0.151 


O.OSO 


Skewness 


0.581 


0.672 


0.687 


0.434 


Kurtosis 


0.304 


1 .040 


1 .055 


0.554 


MiniiDum 


0.94 2 


0.814 


0. S45 


0.723 


Maxiaua 


1.843 


1 .598 


2. C€3 


1.328 



6. Generated 5 


Cth Quantile 


1000 


Replications 




. 


MAP A 


MAP B 


MAP C 


MAP D 


Mean 


24.50 


40.62 


27. Cl 


29. 14 


Standard Deviation 


4.929 


6.758 


5.950 


4.069 


Skewness 


0.54 4 


0.296 


0.449 


0.259 


Kurtosis 


1 .010 


0.233 


C.22C 


0. 1C2 


Minitaum 


1 1.9b 


19.01 


10.51 


18. C9 


Maximum 


50.3 1 


68.74 


49.99 


42. 18 



1000 
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7 



7. 


Generated 


75th Quanti 


le^ 1 000 


Replications 








MAP A 


MAP B 


MAP C 


MAP B 


Mean 




62.7 8 


87.77 


1 H .26 


56. 11 


Standard 


Ce viaticn 


10.37 


12. 12 


12.76 


6.752 


Skewness 




0.41 3 


0.188 


0.244 


0.2C7 


Kurtosis 




0.23 0 


-0.16 


-0. 15 


0.023 


minimum 




36.81 


52. 85 


38.54 


35. 16 


daxioium 




1 10S7 


126.9 


1 13. S 


80.57 


8. 


Generated 


90th Quanti 




Replications 






• 


MAP A 


MAP B 


MAP C 


MAP D 


Mean 




1 15. 6 


145.9 


14 1 .0 


87. 18 


Standard 


Deviation 


16.53 


21.03 


24.61 


11.22 


Skewness 




0.30 3 


0.429 


0.425 


0.326 


Kurtosis 




0.07 1 


0.335 


0.463 


0.369 


tsinimuiii 




68.04 


90.24 


80.58 


54.44 


Maximum 




183. 6 


234.3 


246.6 


126.1 



74 



9 





5§S31i 1 QOC 


Replicat 


ion§ 




MAP A 


MAP B 


MAP C 


Mean 


3 22. 2 


360.4 


397.6 


Standard Deviation 


92.13 


96.13 


117.6 


Skewness 


1.28 3 


1.278 


1.316 


Kurtosis 


2.910 


2.937 


2.728 


Miniaum 


153. 5 


185.9 


190.2 



MAE D 
194.2 
44.55 
0.8S5 
1 . 047 
95. C4 



MaximuiD 



798.5 



889.9 



105S.1 



396.9 
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